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without a singularity. 
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The characteristics of the relativistic cosmology are the expansion and the 
big-bang singular origin of the Universe. The imprint of the latter is believed 
to be seen in the cosmic microwave background radiation observations [1,2]. The 
theoretical model imbibing these features is an exact solution of Einstein's equations 
known as the Friedman-Lemaitre-Robertson- Walker (FLRW) model. It represents 
a homogeneous and isotropic Universe filled with perfect fluid. This is the generally- 
accepted current model of the Universe. Despite its success, there are some questions 
of principles that deserved to be addressed to. The model is highly specialised for 
it is homogeneous and isotropic. These properties could by no means be considered 
generic enough features of the Universe. The Universe should in fact have more 
general initial conditions. Secondly, inhomogeneity may be quite appropriate for 
evolution of large scale structures in the Universe. Above all even at the present e! 
poch, the question of homogeneity 

It is therefore important to find solutions of Einstein's equations when 
homogeneity and isotropy could be dropped. The first step in this direction came 
with the Bianchi models that are homogeneous but anisotropic. For a long time 
only homogeneous models were studied. Inhomogeneity came in through what 
is known as orthogonally transitive G2 cosmologies [3-5]. These are the simplest 
inhomogeneous models that admit two spacelike Killing vectors which are mutually 
as well as hypersurface orthogonal. The class of solutions of Einstein's equations 
of inhomogeneous family was first obtained by Wainwright and Goode [6] and 
subsequently by others [7-9] . Amongst them was a remarkable solution obtained by 
Senovilla [8] which was free of the big-bang (or of any other kind) singularity. It 
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satisfied all the physical reasonableness conditions and had an acceptable equation 
of state p = 3p > 0. Until then on the strength of the singularity theorems [10] it was 
generally believed that the occurrence of big-bang is an inevitable feature of general 
relativity (GR). The conflict was soon resolved [11] when it was shown that the 
solution in question does not obey the assumption of existence of compact trapped 
surfaces and hence the singularity theorems become inapplicable. All prior attempts 
to construct singularity free models had either to ascribe physically unacceptable 
behaviour for matter leading to violation of energy and causality conditions or to 
invoke quantum effects or modiflcation of GR [12,13]. Senovilla's [8] was the first 
singularity free solution, true to GR, conforming to energy and causality conditions. 

The most interesting feature of inhomogeneous cosmology is the non-occurrence 
of singularity in the models. Hence it has no beginning and no end. This is 
aesthetically very attractive and appealing feature, which was first propounded 
by the steady state theory. Ruiz and Senovilla [14] have identified a large family 
of singularity free cosmological models. This family is unique for cylindrically 
symmetric metric separable in space and time in the comoving coordinates (i.e. 
fiuid lines are orthogonal to t = const, surface). In this note we wish to establish 
the uniqueness of the non-singular family by dropping the assumption of cylindrical 
symmetry and in the process we succeed in characterising the perfect fluid models. 

Before we go any further let us note a general result arising out of the following 
two relations [15], 



9,a = - 
2 



(1) 
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where 9, u, a;, iia are the kinematic parameters; expansion, shear, rotation and 
acceleration, Uq, = Ua-iu'^- We have however assumed Ui = ^Jgoo^^- 

We infer from the above relations: 

Lemma : In the absence of shear and vorticity, the expansion of fluid is constant 
over the 3-space orthogonal to the fluid congruence and further, the acceleration also 
vanishes when the quantity g/goo is a separable function of space and time in the 
comoving coordinates. 

Corollary : For the vorticity free spacetime with the separability ( as is the case for 
the metric (3)), acceleration can be non-zero only if shear is non-zero. 

According to the Raychaudhuri equation [16], in the absence of vorticity 
acceleration is necessary for halting the collapse to avoid singularity which in our 
case can only exist if shear is non-zero. Thus non-singular solutions represented by 
the metric (3) will always have to be both inhomogeneous and anisotropic. 

In cosmology in general and G2 cosmologies in particular the metric is generally 
taken to be orthogonal and the comoving coordinates are employed. That is the fluid 
velocity vector is vorticity free and is orthogonal tot = const, hypersurface. Further 
as is common for most cosmological models, we take the metric to be separable in 
space and time coordinates. Hence we write the metric in the form 

ds'^ = Ddt'^ - Adxl - Bdxl - Cdxj (3) 
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where by sepaar ability we mean A = A{t)A{xa) and so on. The velocity field of the 
fluid is Ui = /D5°. 

It can be easily seen that the invariant characterisation of this metric is given 
by (i) 0,a = Oua and (ii) (cr/O), a = 0, i.e. the anisotropy parameter is constant on 
the 3-hypersurface orthogonal to the fluid flow [17]. The condition (i) follows from 
eqn (2) while (ii) can be verified easily by writing a/B for the metric (3). 

Let us couch the general result as a theorem as follows : 

Theorem : Let the metric (3) be a perfect fluid cosmological solution of Einstein's 
equations then it can only be 

(a) if homogeneous, (homogeneity) - Bianchi I or Gg (both homogeneity and 
isotropy) - FLRW model 

(b) if inhomogeneous, G2 (admitting only two spacelike Killing vectors) models 
with or without the big-bang singularity. 

Further the singularity free family as already identified in [14] is unique and is 
cylindrically symmetric. 

Proof : The theorem characterises all the perfect fluid models that the metric (3) 
can represent. But for the separability and orthogonality of the metric we make no 
assumptions. The fluid conditions alone will impose symmetries on the metric. 

The perfect fluid distribution will imply the conditions; Toq, = 0, Tq,^ — for 

To go any further we need the explicit expressions for [18] which look quite 
frightening and formidable. Fortunately, we have discovered an underlying order 
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in them that allows us to write the rest of them from the given two (one each of 
diagonal and off diagonal) by prescribing the appropriate permutation rules. We 
begin by 
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where a subscript denotes partial differentiation and here the assumption of 
separability is not effected. 

The successive cyclic permutations A B C A and 1 — > 2 — > 3 — > 1 will 
give T^,T^ from ; Ti,T^ from T^; and T^,Ti from T/. To write from 
T^, let 0^i2 (i.e. Aq iA2, iT^) and B ^ C ^ D ^ B while follows 

from for 2 ^ 3 ^ 1 ^ ^0 ^ -2{Tl ^ T^) a^nd A ^ D ^ B ^ C ^ A. Thus 
we can write all ten T^, given the two, one each of diagonal and off diagonal. 

Let us begin with the general case where the metric is a function of all the 
coordinates. Eqns. T^o = 0, of which Tiq — reads as 
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D \Bo/B + Co/Cj C ^ \Bo/B + Co/Cj B 
and similarly D2/D and D^/D. It is clear if Aq/A — Bq/B — Cq/C which implies 

(7 = and ■Uq, = (because = ), then the spacetime is both isotropic and 

homogeneous and it can be no other than the big-bang singular FLRW [19]. 

If Aq/A ^ Bq/B ^ Cq/C and since the metric is assumed to be separable in t 
and Xoii the above equation will imply 

Cq/C - Aq/A = hiBQ/B - Aq/A) + miBQ/B + Cq/C) (7) 

and 

Bi/B = kiCi/C (8) 

where ki and ni are constants, and similar equations will follow from D2/D and 
D3/D. We can now set the exact differential for the space dependence of D, 

d{lnD) = {lnD)idxi + {lnD)2dx2 + {lnD)^dxz 

which can be integrated along any path to give the same result. Evaluating it 
along two different paths, we obtain B — C^'^^D — C^^ and A — C^^ . Thus 
the space dependence of the metric is all determined but for the single function 
C{xoi)- Eqns. T^ip = Q iov a ^ (5 determine C{xa) = constant. Hence the metric 
can only represent homogeneous and big-bang singular Bianchi I model. When 
Aq/A = Bq/B ^ Co/C, eqn (6) and its permutants will imply that either it is 
Bianchi I or it admits a spacelike Killing vector. This is the case we consider next. 
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We shall now consider the case of metric (3) admitting a spacelike Killing 
vector, say That means the metric depends only on the two space variables xi 
and X2- 

From (5) it is clear that = tI = T3 will imply the two equations of the 

type, 

where /i, /2, /s are functions of space variables Xa and containing derivatives with 
respect to xi,X2 and X3 respectively. 

In this case /a = in (9). Note that T30 = now and T^o = will lead to 
A{t) = B{t) 7^ C{t) and D{xa) — C^{xa), k= const. Hence a can be non-zero to 
give rise to acceleration which in turn can lead to a viable non-singular case. There 
will be three equations (T12 = plus the two following from (9)) determining the 
space dependence of the metric. 

Since A{t) = B{t), it is possible to perform a coordinate transformation to set 
A.{xa) — B{xa)i which means A = B. K very detailed and involved Lie analysis of 
the equations [20] leads to the conclusion that space dependence can only arise as 
a function of xi -\- X2 or x^ + a;| or Xi/x2 ( i.e. A{xa) = A{xi + X2) etc). In either 
of the first two cases, it could be reduced to the case of one-coordinate dependence 
by coordinate transformations. The last case is obviously singular and would not 
lead to a viable model. 

When Aq/A ^ Bq/B ^ Cq/C, following the above route, eqns. Tao — determine 
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all others in terms of C{xa),a = l,2;Ti2 = fixes C{xa) = {f{xi) + g{x2))^ and 
then eqns. (9) force C{xa) to be a constant. We are again led to Bianchi I model. 

Finally we come to the case of G2 models which admit two spacelike Killing 
vectors that are mutually as well as hypersurface orthogonal. The metric depends 
upon only one space variable. 

Now the spacetime is general enough to sustain inhomogeneous perfect fluid. 
There will occur two kinds of inhomogeneous fluid models, one with singularity and 
the other without it [14] . The singularity free family possesses cylindrical symmetry 
and is unique. That means that the already identified singularity free family [14] 
is unique not only for cylindrical symmetry but for the general metric (3). In this 
case homogeneous models can however occur but will be Bianchi I only. 

Thus is proven the theorem. 

We give below the general enough [21] though not the most general [14] 
inhomogeneous non-singular cosmological models described by 



with a + P = 1. It admits the only two cases, (i) b = c and {ii)b + c = 1 for 
perfect fluid distribution. The former includes the Senovilla's solution [8] when 
c = b = —1/3, a = 1, a = 2, P = —1 and k = 3m; and 



ds^ = cosh^"(A;t) cosh^" (mr)((it^ — dr'^) — cos 



(10) 



m ^ sinh^(mr) cosh^'^(kt) cosh'^'^ {mr)d(p'^ 



15k^ 



cosh-^(A;t)cosh"^(/cr/3) > 0. 



(11) 



p = 3p = 



The latter always gives the stifi' fluid. 
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87rp = 87rp = (6^ - 4)m^ cosh.-^°'-^ {2mt) cosh-^" (mr) (12) 

with a = 1 — 6/2, a = 6(6 — 1), c = 1 — 6, A; = 2m. Clearly 6^ > 4 which means 6 
must lie outside the interval — 2 < 6 < 2. 

Though p > for both b > 2 and b < —2, the geodesic completeness [22] 
demands 6 < 0, and hence the latter case is the truly singularity free solution. At 
the two ends of the range for 6 = ±2, we have the two distinct empty space solutions 
as matter free limits of the stiff fluid. Here again, b = —2 will be geodesically 
complete and it can be thought of as representing field of plane gravitational waves. 
We have also shown [23] elsewhere that a natural inhomogenisation of FLRW open 
model leads to the non-singular metric (9). 

The metric (3) can only describe fluid models of the following kinds; 
homogenous FLRW and Bianchi I and G2 inhomogenous with or without singularity. 
That is inhomogeneity cannot be sustained by any symmetry higher than G2 
[20]. This is an important conclusion; at whatever scale we wish to incorporate 
inhomogeneity in fluid models represented by an exact solution of Einstein's 
equations, the metric has to be G2- Howsoever much it is not supported by 
observations, it appears essential for bringing in inhomogeneity, which cannot be 
ignored even at the present epoch in view of lumpiness in the Universe. The ideal 
solution to the question is if inhomogeneity and anisotropy were scale and time 
dependent. Then at larger scales and at late times the model could have evolved 
to homogeneous and isotropic FLRW. This unfortunately cannot happen for the 
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metric (3) for a/d turns out to be constant for inhomogeneous fluid models. 

It may be noted that the metric (3) is general enough for cosmology. For both 
the assumptions of orthogonality and separability are quite common and shared by 
most of the cosmological models. It is in this context very important that it can 
only admit inhomogeneous perfect fluid models that could never isotropise. The 
only hope is to give up the separability which unfortunately will make the problem 
mathematically formidable. If one can succeed in flnding an inhomogeneous non- 
singular solution of Einstein's equations that isotropises to approximate to FLRW 
as time progresses, it will have very important bearing on our overall cosmological 
perception of the Universe. 
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